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We discuss by analytic means the theory of the high-density limit of the unpolarized two- 
dimensional electron liquid in the presence of Rashba or Dresselhaus spin-orbit coupling. A gen- 
eralization of the ring diagram expansion is performed. We find that in this regime the spin-orbit 
coupling leads to small changes of the exchange and correlation energy contributions, while modi- 
fying also, via repopulation of the momentum states, the non-interacting energy. As a result, the 
leading corrections to the chirality and total energy of the system stem from the Hartree-Fock con- 
tributions. Fhe final results are found to be vanishing to lowest order in the spin-orbit coupling, in 
agreement with a general property valid to every order in the electron-electron interaction. We also 
show that recent quantum Monte Carlo data in the presence of Rashba spin-orbit coupling are well 
understood by neglecting corrections to the exchange-correlation energy, even at low density values. 
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I. INTRODUCTION 

The effect of the spin-orbit coupling in semiconduc- 
tor heterostructures, and more specifically the role of 
the Rashba spin-orbit interaction induced by the asym- 
metry of the transverse confining potential in electronic 
two-dimensional systems^ has attracted in recent years 
great interest. While a main motivation lies in the po- 
tential of new applications,— ~— based on the control of the 
single-particle spin-dependent dynamics through electri- 
cal gating^ - — the corresponding many-body problem is 
of fundamental relevance and not yet fully investigated. 

The two-dimensional electron liquid, in the presence of 
Coulomb interaction and a rigid neutralizing background, 
is a classic problem in solid-state physics, if the simplest 
effective mass approximation is assumed for the kinetic 
term»ii On the other side, concomitant band-structure 
effects have often large observable consequences ! 10 ' 12 : 13 
The special form of a generalized spin-orbit coupling ap 
plicable in a number of cases is described in Ref. 
together with a detailed analysis of the exchange en- 
ergy. Several other aspects of the electron-electron in- 
teraction in the presence of spin-orbit couplings were 
addressed in Refs. Il5l - [32l For example the quasipar- 
ticle propertieS ) 15 : 16 ^ 3 -^"— ^ the Hartree-Fock phase- 
diagram ; 17 : 22 ' 23 the spin susceptibility ) 19 : 23 ' 30 ' 32 and the 
plasmon modesi 18 ' 21 ' 24 ' 31 ' 32 We restrict ourselves here 
to pure Rashba^ or equivalently Dresselhaus ; 33 ' 34 spin- 
orbit coupling. In this case, quantum Monte Carlo data 
for the total energy were recently obtained in Ref. [25|. 
Numerical results for the total energy also appear in 
Ref. |29| . within the random phase approximation (RPA) . 
We focus in the following on the high-density limit, when 
the effect of the Coulomb interaction can in general be 
studied perturbatively. 

The leading correction to the non-interacting energy 
is the exchange contribution, while higher order terms 



correspond to higher powers of the standard density pa- 
rameter r s [see Ref. [TT| and Eq. (O below] . The first two 
terms of the small r s expansion of the correlation en- 
ergy in two-dimensions are well known in the absence of 
spin-orbit coupling^ They consist of the second-order 
correlation energy and a r s In r s contribution which is 
obtained as an infinite sum of diverging ring diagrams. 
The elegant resummation procedure patterns the treat- 
ment of the leading In r s correlation energy in the three- 
dimensional case^ An exact formula for the polarization 
dependence of the r s In r s contribution was recently de- 
rived in Ref. 

Here, an additional dependence on the (dimensionless) 
Rashba coupling a is present. However, the strength of 
the spin-orbit interaction is more appropriately expressed 
in terms of a parameter g which is proportional to r s [see 
Ref. [3 and Eq. (2J) below] . Hence, an additional density 
dependence is introduced by g. This makes the effect 
of the spin-orbit interaction small, since the correction 
to the total exchange-correlation energy is multiplied by 
a factor at most of order g 2 (i.e. an even power in the 
spin-orbit coupling), which is vanishingly small at high 
density. Therefore, an accurate result is obtained per- 
turbatively. In fact, a general argument for the energy 
expansion was derived in Refs. [23| and[28| and implies that 
the g 2 term is actually vanishing (see also Ref. |29I ). 

In this paper, we analyze how the high-density expan- 
sion of the energy is modified in the presence of spin- 
orbit coupling. The explicit analytic form of the lead- 
ing exchange-correlation correction is obtained in the fol- 
lowing and found indeed to be proportional to g 4 lng, 
from the exchange energy. The second-order correla- 
tion energy is studied numerically, as in the case with- 
out spin-orbit interaction, and the extension of the ring- 
diagram sum is also discussed, and shown to display a 
non-analytic behavior in the limit of small r s and g. Cor- 
rections to these higher order contributions are also found 
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to be of higher order than g 2 , as expected, and can there- 
fore be usually neglected. Finally, while the main body 
of the paper is devoted to the asymptotic expansion at 
small r s (a regime often relevant for heterostructures with 
large spin-orbit coupling), we also analyze the quantum 
Monte Carlo results of Ref. ,25, which are all at r s > 1. 
We propose here a simple interpolation formula for the 
energy which is in remarkable agreement with the numer- 
ical data. 

The detailed outline of the paper is as follows: in Sec. 
[IT] we formulate the problem and establish our notation. 
We review the properties of the non-interacting system 
and the known results for the exchange-correlation en- 
ergy in the absence of spin-orbit coupling. We also define 
here the corrections to the exchange-correlation energy 
of the electron liquid due to the spin-orbit interactions, 
which are the main focus of our work. We first show in 
Sec. IIIII that such corrections are generally small, by re- 
minding the reader about some useful exact properties of 
the perturbative expansion; 23 ! 28 and by an analysis of the 
quantum Monte Carlo results of Ref. [25|. We then explic- 
itly determine such corrections in the high-density limit 
of the electron gas. The exchange energy, the second- 
order correlation terms, and the classic ring expansion of 
Ref. [H are revisited and extended in Sec. [TV] where we 
obtain the change of the exchange-correlation energy and 
of the momentum space occupation to leading order in 
the spin-orbit coupling. Both analytical and numerical 
results are provided, which are summarized in Sec. [V] 
In this last section, an alternate physical limit is also dis- 
cussed. Finally, the details of a number of calculations 
have been provided in Appendices O [B] and [C] 



II. 



FORMULATION OF THE PROBLEM 



The system is described by the hamiltonian 



* = E^ + jE 



(1) 



where terms related to the presence of a uniform neutral- 
izing background have been omitted for simplicity. The 
single particle operator Hq is given by: 



H 



2m 



+ a (a x p y - 6yp x ), 



(2) 



where we consider electrons confined in the (a;, y) plane 
and cr x (y) are Pauli matrices. The spin-orbit term is usu- 
ally referred to as a linear Rashba spin-orbit coupling and 
is generally present when the confining potential in the 
z direction of a quantum well is asymmetric. 1 An equiv- 
alent term, the Dresselhaus spin-orbit coupling, arises 
instead for a lack of inversion symmetry in the crystal 
structured We consider Eq. (|5J) as a model case, al- 
though a similar analysis can be carried out for other 
types of spin-orbit interaction, relevant in other experi- 
mental cases i 10 ' 12 i 14 



We use in the following dimensionless units. The prop- 
erties of the electron liquid are completely determined by 



a 



ha 



and 



1 



(3) 



where a is a dimensionless spin orbit coupling and r s 
the usual density parameter, with as the effective Bohr 
radius and n the number density. It is also useful to 
introduce the following dimensionless coefficient g 



= V2c 



(4) 



which better than a represents the strength of the spin- 
orbit term. In fact, g is approximately equal to the ratio 
of the spin-orbit energy to the kinetic energy, which are 
respectively proportional to g/r 2 and 1/r 2 [see also Eq. 
(|10|)]. The wave vectors are expressed in terms of the 
Fermi wave vector kp = \j2im and the energies are in 
Ry units. Finally, throughout the paper we often use the 
notation p = k + q and p' = k' — q (so that q will not 
explicitly appear in many expressions). 



A. Non-interacting electrons 

The non-interacting problem is completely determined 
by g. The eigenstates of Hq can be written as 



¥>k,±(r) 



ik-r 



(5) 



where L is the linear size of the system and (/>k is the 
angle formed by k with the x axis. The eigenstates have 
spin quantized perpendicular to the wave vector k, as de- 
scribed by the spinor functions |k±) [which are defined by 
Eq. (|5|)]. The corresponding eigenenergies (in Ry units) 
are equal to \ e±(&) where 



e±(fc) = k 2 Tgk. 



(6) 



We also define the generalized chirality which de- 
termines the occupation functions n±(k) in momentum 
space i 17 ' 22 The two relevant regimes are depicted in Fig. 
Q] At high density (second panel of Fig. [T]), two bands 
are occupied and x coincides with the regular chirality: 



X 



N+-N- 
N + +N- 



, where N± is the total number of electrons 



in each band. The occupation takes the form 



n ± (k) = e(y/l±x-k) forO< X <l, (7) 

where 9{x) is the usual step function. At low density 
(first panel of Fig. Q]), the higher band is empty and the 
occupation function reads 



(k) = 0(^x~TT-k)-9(^ x ~T-k) 1 (8) 



where \ ^ !• We note that Eq. ([8j corresponds to a ring 
in momentum space and the regular chirality is 1 in this 
case, irrespectively of the precise form of the occupation. 
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FIG. 1. Two different ways of occupying non-interacting chi- 
ral states in k space. 



The Fermi surfaces are determined in all cases by the 
radii 



(9) 



The non-interacting energy at generic \ is expressed 
as follows (when \ < 1) 



£o(g,r s ,x) = — 5 2 o ■ ( 10 ) 



The first term is the kinetic energy, and one has to replace 
1 + X 2 with 2\ if x > 1. The second term is the spin-orbit 
energy. The non-interacting ground state is specified by 
the value of x which minimizes Eq. (jlOp for given values 
of a and r s and is therefore uniquely determined by the 
parameter g of Eq. (Q} 



Xo(fiO 
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1 - £ for < ff < 72 
+ £ for .g > V2. 



(11) 



The corresponding ground state energy is obtained ac- 
cordingly: 



£o(g,r s ) = £o{g,r s ,xa(g)) 



B. Exchange-correlation energy 



(12) 



The exchange-correlation energy of the electron liq- 
uid without spin-orbit coupling is a relatively well known 
quantity^ The perturbative expression at high density 
reads 3 ^ (for the unpolarized case) 



£xc( L Ts) 



Vv 2 0.385-^(10-3^)r s lnr s 



371-7% 



3tt 



where the first term is the exchange energy. The constant 
results from the numerical integration of the second order 
correlation energy, and the last contribution is obtained 
from the infinite sum of diverging ring diagrams, simi- 
lar to the original calculation for the three-dimensional 
case^ At generic values of the density, £ xc (r s ) is ob- 
tained numerically with the Monte Carlo metho d 11 ' 38 
(see however Ref. [37| for the polarized case). 

On the other hand, the exchange-correlation correction 
due to the spin-orbit coupling is to date not accurately 
known. We introduce the following definition 



£(g,r s ) = £ (g,r s ) + £ xc (r s ) + 5£ xc {g,r s ) 



(14) 



where the non-interacting energy is given by Eq. (|12[) 
and S£ xc (g 1 r s ) represents the correction to the exchange- 
correlation energy associated with the spin-orbit cou- 
pling. The latter is generally neglected, for exam- 
ple in density functional studies including spin-orbit 
interactions ! 39 ' 40 A partial justification to this procedure 
is given in Ref. 28, which shows that this correction is ac- 
tually vanishing to quadratic order in g for the particular 
case of the Rashba or Dresselhaus spin-orbit interaction. 
However, such a correction is not zero in general and can 
reasonably lead to important effects in the case of large 
g (e.g., at low density) or for other types of spin-orbit 
interactioni 1 ^ 

Formally, under the assumption that the system be- 
haves as a Fermi liquid, the total energy of the inter- 
acting system can be obtained as a perturbative expan- 
sion (see next section) constructed from a particular non- 
interacting state, as for example the one used to obtain 
Eq. (fTU]) . Therefore, in the case with spin-orbit coupling, 
the total energy acquires an additional dependence from 
the chirality x °f the non-interacting state used in the 
perturbative expansion. This does not need to be the 
starting non-interacting ground state. As a consequence, 
we can quite generally write the total energy as 



£(g,r s ,X) = £o(g,r s ,x)+£xc(r s ) + S£ xc (g,r s ,x) , (15) 



(13) 



where £o(g,r s ,x) is given by (flU]) . The dependence of 
the total energy on x, a t given values of g and r s , is also 
obtained in the Monte Carlo study of Ref. H^, where x 
corresponds to the occupation of the initial trial wave 
function. The data are reproduced in Fig. [21 

The actual value of the interacting generalized chi- 
rality x{g> r s) is obtained by minimization of Eq. (|15p . 
which also yields the corresponding ground state energy 
Eq. (fT?]) . It is important to realize that there are the 
two different contributions to S£ xc (g,r s ). The first one 
is given directly by S£ xc (g, r s , x(g, i" s )), while the second 
one arises from the renormalization of x (i- e -, the repop- 
ulation) in the non- interacting energy £o{g, f s ,x)- 
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III. FORMAL PROPERTIES OF THE 
DIAGRAMMATIC EXPANSION 

We report here, in view of their usefulness, two exact 
results concerning the perturbative expansion of the en- 
ergy and the quasiparticle self-energy. These results have 
been obtained in Refs. andHH for a generic two-body 
potential v(q). The first of the two results pertains to all 
diagrams D contributing to the total energy and reads: 



d 2 D 



d 2 x 



d 2 D 



d 2 g 



d 2 D 



dxdg 



(16) 



This allows one to infer that, for small g and Xi the total 
correction to the exchange and correlation energy must 
behave like 



SS xc (g,x)=C(g- X ) 2 + -- 



(17) 



where in general C is an unknown constant. 

The second result concerns the self-energy E M (fc,cj) 
which is seen to satisfy a similar exact relation to lin- 
ear order in g 



dg 



H dE (fc,u) 
2 dk 



(18) 



where So(fe, u>) is the (interacting) self energy in the ab- 
sence of spin-orbit coupling (i.e. g — x — 0). This rela- 
tionship allows us to write: 



£u(fe,u>) = S (fc - 



MS 
2 ' 



w) + .. 



(19) 



1 ± | are the same as in the 



From this formula we conclude that to linear order in g 
all quasiparticles properties (e.g., the lifetim o 16 i 26 ' 27 ) on 
the Fermi surfaces k± 
absence of spin-orbit coupling 

For applications of these exact results we refer to 
Ref. [HI, which shows the validity of Eq. (fH?)) within the 
RPA, and to the numerical calculations of the quasipar- 
ticle lifetime of Refs. [26| and[27]. The total energy was ob- 
tained in the quantum Monte Carlo study of Ref. [25| and 
we present in Fig. [5] their numerical data, together with 
the simple approximation of setting 8£ xc {g, r s , x) = in 
Eq. (p~5|) . The value of the exchange-correlation energy 
£ xr (r s ) m the absence of spin-orbit coupling is taken from 
Ref. [33- As seen, not only are the data at higher density 
(r s = 1) in very good agreement with the curves obtained 
neglecting 8£ xc {g, r a ,x), but also the low-density data are 
well described by this approximation. This implies that 
the spin-orbit coupling shift of the exchange-correlation 
energy is generally small. Notice that the accuracy of the 
numerical data does not allow us to extract the constant 
C of Eq. (fTT|) as function of r s . This we were able to 
obtain analytically in the limit of small r s , as described 
in the following sections. 
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FIG. 2. Numerical data (solid dots) from Ref. [25| for the 
total energy £(g,r s ,x) of Eq. (fl5|) . as functions of x and at 
different values of g and r„. We obtain the solid lines by 
setting S£ xc (g,rs,x) = in Eq. (|15|) and using the value of 
£xc(r s ) of Ref. |38|. For reference, the empty dot in the top 
panel is the noninteracting energy without spin-orbit coupling 
(at g = x = and r s = 1). 
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FIG. 3. Correction to the exchange energy (in Ry units at 
r a = 1) for the unpolarized state with Rashba spin-orbit cou- 
pling, at finite values of \- To obtain the value for generic r a 
one simply divides by r s [see Eq. (|20[l ]. 



form of the exchange correction for the specific case of 
Rashba or Dresselhaus spin-orbit coupling. This is ex- 
pressed as follows^ 



(20) 



where 8£ x (x) is plotted in Fig. [3J The expansion of 
S£ x (x) at small x is given by 



if,_X , 23 



6£ ^ = W^ ln f + l2 



(21) 



where higher order terms are C(x 6 )- The term quadratic 
in x is missing, in agreement with Eq. (|17p (with C = 0). 

By making use of Eqs. flQ"]l. (|T3)>. (|2U|) . and dHJ), one 
obtains that the energy is minimized when 



IV. PERTURB ATIVE CONTRIBUTIONS OF 
THE COULOMB INTERACTION 

We examine in this section how the first terms of the 
high-density expansion are modified by the spin-orbit 
coupling. In particular, we discuss the exchange energy, 
the second-order correlation energy, and the sum of the 
ring diagrams. One has to notice that, if the bare value 
of the spin-orbit coupling a is kept constant, the r s — > 
limit also corresponds to a vanishing strength of the spin 
orbit coupling g. This is clear from Eq. (jj) and is sim- 
ply understood as follows: at high density the spin-orbit 
energy grows like ahkp, but becomes negligible with re- 
spect to the kinetic energy, which is proportional to k 2 F . 

Furthermore, since x(g, r s ) is given in first approxima- 
tion by the non-interacting expression Eq. this limit 
corresponds also to a vanishing value of x — 9- There- 
fore, we will obtain an expansion of Eq. (|15[) in the small 
parameters r Sl g and Xi which are all of 0{r s ). 

Using the result of the previous section at small g 
and x one can infer that to a generic contribution of 
order C(r") (in the absence of spin-orbit coupling) cor- 
responds a correction cx r"(g — x) 2 which is vanishing to 
lowest order for the ground-state energy. Therefore, the 
leading analytic contribution to 5£ xc (g,r s ) in Eq. (|14[) is 
C( r ™# 4 ) = O (r" +4 a 4 ). This argument is valid for the 
exchange energy and the second-order correlation terms. 
On the other hand, due to their non-perturbative resum- 
mation of all orders, the series of diverging ring diagrams 
requires a more careful analysis. 



A. Exchange energy 



x(g,r s ) = xo(g) 



24tt 



r s g 



8 6 



(22) 

an expression that represents the analytic form of the 
small enhancement of the chirality numerically obtained 
in Ref. [Til The relative correction is of order 0(rf) and, 
as anticipated, it is quadratic in the spin-orbit coupling. 

In calculating the correction to the total energy, one 
has to notice that Xo{g) is a stationary value of the non- 
interacting energy, and therefore the corrections to the 
non-interacting energy due to the renormalized value of 

X are of order O ^ 9 r ^ ^ = C(rf ) and can be neglected. 
The leading term is therefore given by 



S£i x c Hg,r s ) = 



487rr s 



8 12 



(23) 



which represents the leading contribution to the 
exchange-correlation energy correction. We discuss next 
the higher order correlation terms. 



B. Second-order correlation energy 

The second-order correlation energy £2(5, x) is ob- 
tained by standard perturbation theory. In the inter- 
mediate state two electron- hole pairs are present, such 
that occupied states with wave vectors k, k' and chiral 
indexes fi, p! have scattered to new unoccupied states 

(k )A i)->(p,i/) and (kV)^(p>'), (24) 



The exchange energy is the main contribution of the 
electron interaction at high density. The Hartree-Fock 
approximation of the two-dimensional electron liquid in 
the presence of spin-orbit coupling was already studied 
in our previous worki 14 i 22 i 23 We derive here the explicit 



where p = k + q and p' = k' q. As it is well known, 
there are two different ways to scatter back to the original 
states. For direct processes 



(p, v) -> (k, fx) and (p', v') -> (k', //), 



(25) 



6 



which gives 



- r 

1-3 



rfq 



47T 3 4^ ../ g 2 



dk / tflc' 



x n M (fc)n M ,(fcQ(l - n„(p))(l - rv(pQ) 
e w (p)-e M (*0+<v(p')- V( fc 

x|(pHMi 2 |(pV|ky>| 2 . (26) 

For exchange processes 

(p,iO->(kV) and (p',^')-^ (k,/i), (27) 



which corresponds to 



dq 



q\k - k' 



n^k)n^(k')(l - n„(p))(l - rvfc/)) 



dk / dk' 



e„(p) - e M (fc) + e„/(p') 



(28) 



x(pi/|k M ) (k/ilpV) (pV'|kV) (kVlpt 
Finally, the total second-order correlation energy is 

£ 2 (5,x) = £ 2 D (.9,x) + £ 2 x (s,x) ■ (29) 

As in the case without spin-orbit coupling, Eqs. (|26p 
and (|28p cannot be evaluated analytically in general. Fur- 
thermore, computing these multi-dimensional integrals 
with spin-orbit coupling is complicated by the presence of 
singularities in the integration domain, from the energy 
denominators. In fact, the excitation energy is guaran- 
teed to be positive when x — XoG?) but in the general 
case the energy denominator can be zero or negative. By 
restricting ourselves here to the simplest case x — Xo(<?)j 
which is correct to leading order at high density, we plot 
in Fig. 2] the direct and exchange second-order integrals 
as functions of g for g < y/2 (corresponding to x < !)■ 
The sum of the two is plotted in a wider range of values 
in Fig. [5j We notice in Fig. [4] that in the limit g both 
functions £® and £^ display a flat behavior in agreement 
with the vanishing of the the g 2 contribution. The direct 
term is larger and dominates the sum as displayed in the 
inset of Fig. [5] The characteristic behavior, similar to the 
case of the exchange energy (see Fig. |3J), suggests a g 4 In g 
leading term. It is also remarkable that at large values 
of g the correlation energy diverges. This limit of large 
spin-orbit coupling is highly non-perturbative, as already 
revealed by the Hartree-Fock treatment. Within that ap- 
proximation, the non-interacting states are strongly dis- 
torted by the Coulomb interaction and form special spin- 
textures in momentum spacei i 22 ' 23 

Analytic formulas for the second-order correlation en- 
ergy can be obtained at small g and x- Expanding to 
second order yields 

(g-x) 2 , 



£?{g,x) -- 0-614 - 



0.229 



(30) 

(31) 
and C 



i n agreement with Eq. (jTT)) (with C = 
respectively). The explicit calculation is detailed in 
Appendix [X] 
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FIG. 4. Plot of the second-order correlation energy for the 
non-interacting ground state, as function of g. The general- 
ized chirality is \ = Xo(ff) and the range of both plots is such 
that g < v2, which gives x < !• The top panel shows the 
direct term £F (g,Xo(g)) and the lower panel the exchange 
term (g,Xo(g))- The points represent numerical results 



from Monte-Carlo integrations of Eqs. 
solid lines serve as a guide for the eye. 



26j and J2HJ, and the 




g/V2 

FIG. 5. Plot of the total second-order correlation energy 
&(ff, Xo(<?)), with non-interacting occupation, as a function 
of g. The inset shows the region of small g, when x < 1> an d 
is obtained as the sum of the curves displayed in the previous 
Fig. |4] The points represent numerical results, and the solid 
lines serve as a guide for the eye. 
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C. Ring diagrams 

The higher order terms (n > 3) in the perturba- 
tive treatment of the two-dimensional electron liquid are 
in general diverging for the bare Coulomb interaction. 
However, a method to obtain the next leading correc- 
tion to the correlation energy was devised for the three- 
dimensional case.— It consists in summing to infinite or- 
der the (regularized) most diverging diagrams so that 
the final result is finite^ This method was applied in the 
two-dimensional case in Ref . |35| and is extended here with 
suitable modifications to include the Rashba or Dressel- 
haus spin-orbit coupling. 

The expression of the ring diagrams reads (n > 1 ) 



+oo />oo 

du I q l dq . _ , 
I J-oo Jo V 2\ / 2nq J 

(32) 



where Q q {u) is given by 



^2 [ dk n M (fc)(l-n v (p))^ 



( £ „( P )- £(1 (fc))|( P HMp 



(p)-e M (fc)) 2 /4 + U V 



(33) 

£^ (g, r s , x) uas a (formal) dependence on r™~ 2 . Except 
for the n = 2 term, which is merely a compact formula 
for Eq. (j2"6"]l . these expressions diverge. Summing them 
up to infinite order we arrive to 



1 



du 



In 1 



Q q (u)r s 
2V2ti q 



q dq 



1 ( Q q (u)r s Y 

2 V 2V2irq J 



Q q {u)r s 



2V2nq 



(34) 



Rather than plunge into a numerical analysis, we en- 
deavor next to extract the analytic behavior of this con- 
tribution on the variables r s , g and \- We begin by as- 
sessing the behavior of the function Q q (u) for g, x 0. 
We first notice that the q — > limit of this function is not 
analytic. To see this consider that, as shown in Appendix 
IB1 for fixed q and g,x — > 0, Q q (u) behaves accordingly 
to the general form given by Eq. (fTTj) . i.e., like 



» ~ Q(°>(u) + 



where Q q °\u) is the value of Q q (u) for g = x = and 



(5-X) 2 + 



(35) 



Q 9 q 9 (u) 



9 2 Q q 



dg 2 



is given in Eq. (TBI 



On the other 
-> for fixed 



hand this relationship does not hold as q 
g and %. I n this case one can neglect in Eq. (|33|) terms 
involving scattering to the opposite branch. This is justi- 
fied since the factor |(z/p|/ik)| 2 is ~ 1 for the intra- band 



and 



q for the inter-band contributions. Then 



^ fdk 2n M (fc)(l 



n^p)){k- if) cosfa 



M )2 , 



k + k + 



2vrfc_ 



(36) 



where ii(tt) = 1 
X < 1 so that fc± 



l/yl + I/" 2 and we have assumed 



y/1 ± x- We have also defined k± 
k± =F g/2. In particular, by setting % = g, we obtain 



? (U) =2 47Ti?( U ) 



irug 



2(1 + u 2 ) 3 / 2 



(37) 



which of course is in violation of Eq. dT7|) 4^ 

As a consequence one cannot immediately infer that 
Eq. ([34]) does in fact satisfy the general property Eq. (|17[) . 
This however proves not to be a problem since Eq. (|36[) 
only applies within a region of q values of extension much 
smaller than g and x- Accordingly, in order to obtain the 
leading contributions to £n(g,r s ,x) when x-,9 ~ ► one 
can safely make use of Eq. ([33]) for Q q (u). Therefore, the 
correction to the generic ring diagram can be formally 
written as 



«<») 

27rr„ 



+co />oo 

~ ' / 9 2 e?9 

oo JO 



2\/27rq 



Q£ 0) (u) 



(38) 



to lowest order in g and x- Summing the geometric series 
(for n > 2) we have 



16tt 3 



{g-xf 



(«) Qf(") 



2 v / 27rg + Q^ 0) (u)r s ' 
(39) 

where the integral is is logarithmically divergent at r s — > 
0, due to the small q integration region. Therefore we 
can approximate 



5£ r 



■ In? 



2\/2(27r)' 



-{g-xf 



where Q 



(0), 



(40) 

AttR(u) [with i?(it) defined below 
Eq. (36])] and Qg 9 (u) ~ -3tt|u|/(1 + it 2 ) 5 / 2 which results 
from the small q limit of (|B1[) . We can then integrate 
Eq. (gOJ) to obtain 



8£ii(g,r s ,x) 



56 - 15?r 
40\/27r 



(.9 - X) 2 + 



j In? 



(41) 

where as usual we neglect higher order terms in g and x- 
A discussion of the physically different limit of small 
r s and finite g, x? f° r which Eq. (|36[) becomes applicable, 
is provided in next section. 



V. DISCUSSION 

By adding up the contributions from the previous sec- 
tions we obtain that the total energy per particle in the 
presence of spin-orbit coupling of the Rashba or Dres- 
selhaus type has the following form, in the limit of high 



8 



density and small values of g and x 

= £o{g,r s ,x) 

8^2 , V2 X , 23 



X ( In 

3tt 48tt V 8 12 



0.385 



2y/2 



(10 - 3tt) 



56 - 15tt 
40V^vr 



(5-x) 2 + 



In? 



(42) 

where the non- interacting energy £o(<?j r siX) i s given in 
Eq. (ITOl) . We included in Eq. (1421) all the quadratic terms 
in x an d g as well as, for the exchange energy only (third 
line), the term of order C(x 4 ) which represents in prac- 
tice the leading correction at high densities. The fourth 
line shows the contribution of the direct and exchange 
second order diagrams, from Eqs. (|3"0|) and (|3"Tj) . Finally, 
the last term (in the fourth and fifth lines) represents 
the ring diagrams correlation energy. For g = x = 
this expression recovers the classic result for the two- 
dimensional electron liquid of Rajagopal and Kimball of 
Ref.Hl 

From Eq. (l4"2"j) it becomes clear that, for small r s and 
given spin-orbit coupling a, the largest effect of the in- 
teraction is due to the exchange energy. In particular, 
we could obtain in Eq. (|2"2")l how the best possible re- 
population, i.e., the value of x m equilibrium, is modi- 
fied by the interactions to leading order. We also deter- 
mined in Eq. ([23| the leading correction of the exchange- 
correlation energy due to the spin-orbit coupling. In 
agreement with the general property discussed in Sec. Mil 
such corrections are very small. Beyond the pertur- 
bative regime, we have shown (see Fig. [2]) that recent 
quantum Monte Carlo data are well reproduced by ne- 
glecting the correction to the exchange-correlation en- 
ergy due to the spin-orbit coupling. Our results lend 
some measure of comfort to the otherwise uncontrolled 
procedure of making use of spin-orbit coupling free ex- 
change and correlation functionals within density func- 
tional calculations . 39 ! 40 

We end our discussion with an analysis of the formally 
interesting case of the limit of small r s at finite g and x- 
This corresponds to the case of a diverging bare spin-orbit 
coupling constant a [see Eq. (j4])], a scenario that does 
not apply to the plain two dimensional electron liquid 
treated up to this point. In this case the otherwise small 
'anomalous' integration region contributing to Eq. (|34[) 
becomes dominant and Eq. (|36j) is the appropriate form 
for Q q (u). While the reader is referred to Appendix ICl 
for the details of the calculation, we give here the result 
for the correlation energy to quadratic order in g and x 

d£ R (g,r s ,x) ■ 



,208 51ti\ 9 



A2 9n, 9 27, 



(43) 



This expression is distinct from Eq. (|4Tj) and violates the 
general form cx (x~g) 2 - As a consequence, the quadratic 
term in the spin-orbit coupling survives and the correc- 
tion is proportional to g 2 r s lnr s instead of g i r s lnr s (us- 
ing x — g)- F° r the exchange energy and the second- 
order correlations the non-analyticity in the r Sl g — > 
limit is not present, and the same results obtained before 
are valid here. Despite this fact, the ring-diagrams con- 
tribution remains a subleading correction since Eq. (|43[) 
is applicable only when g 3> r s which clearly implies 
g^hig/rs ^> <7 2 r s lnr s , i.e., the exchange energy correc- 
tion is larger. 

As a final remark, we notice that the nonanalyticity 
of the r s , g — > limit discussed here becomes relevant in 
gated heterostructures since the increase of the density is 
naturally accompanied by a modification of the confining 
potential. If smaller values of r s require higher values of 
the external electric field (e.g., proportional to the elec- 
tron density^) , larger values of the spin-orbit coupling a 
are obtained at the same time, thereby making the limit 
of Eq. (|43p meaningful. Furthermore, in heterostructures 
with a back gate, r s and a can be controlled separately. 
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Appendix A: Derivation of Eqs. (|30p and (|31|) 

We consider here the second derivatives of S-^ig^x) 
and £2{g,x)i given by Eqs. (|26l and (|28|) . respectively. 
We start by listing some useful formulas for spin summa- 
tions: 

N 

n^ pj ' ij i pi + i ^ i + i ) = 2 ' ( ai ) 



N 



5z n^ pi/ii i pi + i ' ii + 1 ^ = 2 " 



pjPk 



(A2) 



from which the following three identities are obtained 

XI ^(vp - mokphmi 2 = 2k q 



k ' 



(A3) 



J2 |<pHkM)l 2 KpV|ky) 



*{vp- iik + v'p' - ii'k'f = 8q 2 . (A4) 

X (pHM (k M |pV}< P V|ky) (kyip*,) 



x{vp - [ik + v' p' - n'k'f = 0, (A5) 
where, as in Eqs. (121)1) and (|2"8"|) . p = k+q and p' = k' — q. 
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We now examine 



After calculating the second 

derivative of the integrand, which only involves the en- 
ergy denominator, the spin summation can be evaluated 
by making use of (|A4I) . One can next evaluate the angu- 
lar integration in the dq integral, and than integrate in 



dk 



dk' y to obtain the following expression 



2cD 



d 2 £ 



dg 2 



\r%f dk X d K L ^ )L ^ 



^ Jo T J-i (<l + kx-k' x ) 3 

(A6) 

Here L(q,k x ) — J no(k)(l — no(p))dk y , which in the in- 
terval — 1 < k x < 1 gives 

!0 if q < 2 and k x < -§, 

£(k x ) iik x >-q+l, 
£(k x ) — £(k x + q) otherwise. 

(A7) 

where we defined £(k x ) = 2-Jl — A; 2 .. We finally have 
evaluated the integral (|A6|) and obtained a result numer- 



ically equal to — ^. The remaining second derivatives of 
£2 can be obtained from the general relation Eq. (fTS)) 
but, as an example, we calculate them here explicitly. 



The expression of 



dgdx 



can be simplified by making 



use of (jAT]l 
we obtain 



d 2 £? 



and (| A2 j> . After some further manipulation 



1 

47T 3 



dq 

I 1 



d\dg 

x [(k.q)(l-no(p))«(l-fc) 



dk / dk' 



n (k')(l-n {p')) 



(<? 2 +q-(k-k')) 2 
(p-q)n (k)6(l-p)]. 

(A8) 



Then, after a change of variable k — > k q = p in the 
second term in the integrand, angular integration in the 
dq integral, and integration in dk y and dk' y , we obtain 



d 2 £? 



dxdg 



1 



% ( dk x d K k ^- k '^ 
q 2 J-i 



no(p) 



no(p) 



(q + k x - k' x ) 2 (k x + k' x f 



(A9) 



we obtain 



d 2 £? 



dx 2 



d 



1 

^ Jo '/ 

1 - n {p) 



dk / dk'n (fc')(l -n (p')) 



dk q 
x 5(1 -p) - 



k x k 
2 cos(</>k 



- fc ) - 



d 



n (k) 



^dp q + k x - k' x 
^ S(l-k)S(l- P ) 



(All) 



where in the dk, dk' integrations, the x axis is chosen 
to be along q so that p = \Jk 2 + q 2 + 2qk x . It is also 
convenient to change variable k — > — k — q in the sec- 
ond term in the square brackets. The contributions from 
the first two terms in the square brackets involving the 

derivative dn g^ = —5(1 — _p)ff cancel exactly the third 
term. In fact, we obtain that the coefficient, multiplying 
the product 5(1 — k)5(l — p) has the following form 



dp/dk 
q + k x - I 



dp/dk 2cos( 



q + k x 



(A12) 



which can be simplified using k x — cos 0k = —q/2. Fur- 
thermore, dp/dk = (k + q cos 4>k)/p = 1 — q 2 /2 and 
cos(0k — </>p) = 1 — q 2 /2 as well. Therefore (|A12[) is seen 
to vanish identically. The term surviving in Eq. (|A11[) . 
after integrations in dk y , dk' y , is the opposite of (|A9|) , 
Therefore, we can summarize the final result as follows: 



d 2 £? 



dg 2 



d 2 £? 



dxdg 



d 2 £? 



dx 2 



(A13) 



We turn now to the second derivatives of £^ ■ The fact 



that 



= is immediately obtained from the spin 



dg 2 

summation (| A5|) and the mixed derivative is also found to 
vanish upon spin summation. Finally, the second deriva- 
tive with respect to x 1S found to be vanishing after rather 
cumbersome manipulations. These results are consistent 
with the general property (TIBI and we conclude that: 



d 2 £? 



dg 2 



d 2 £? 



dxdg 



d 2 £? 



dx 2 



(A14) 



where p = \Jl + q 2 + 2qk x . Finally, Eq. (|A9j) can be 
transformed to the opposite of (|A6|) , by means of an inte- 
gration by parts in dk x of the two terms of the integrand 
and and a suitable change of variable k x — > — k x — q in 
the second one. 



The last term is 



~Sx r 



It is convenient to trans- 



form terms containing double derivatives of the occupa- 
tion functions in the following way 



/(k) 



d'n^k) 



dx 2 



dk 



dk 



5(l-k)dk. (A10) 



After then performing the spin summations explicitly, a 
change of primed and unprimed variables to collect simi- 
lar terms, and the angular integration in the dq integral, 



Appendix B: Second derivatives of the ring 
diagrams contribution 

We consider here the second derivatives with respect 
to g and x 01 Qq( u )i defined in Eq. We begin by 

noting that 



d 2 Q q (u) 



dg 2 



= 16 



dk 

— n (k)(l - n (p))(q + 2k x ) 

q 



(q + 2k x ) 2 - 12u 2 
[(q + 2k x ) 2 + Au 2 ] 3 ' 



(Bl) 



where we have set x = 9 — and performed the spin 
summation (which gives a factor of 2). Consider now the 
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mixed derivative 
d 2 Q q (u) 



d\dg 



/i<5(l - k) 



(l-«o(p)) 



-p) 



n (k) 



(vp - Hk)\(pv\kfi)\< 



4[(q + 2k x ) 2 - 4u 2 ' 
: [(( ? + 2fc c ) 2 +4 W 2 ] 2 



(B2) 



By change of variable k — > — k q and relabeling n v 
in the second term of the integrand (from the large square 
parenthesis), we can cancel the no(p) contribution in the 
first term. Finally, the spin summation can be performed 
by using Eq. (IA3|) and the integration in <ik gives 



d 2 Q q (u) 



dgdx 



4cos0 (g + 2cos( 



4u 2 



q [(q + 2 cos 4>) 2 + 



212 



(B3) 



The last second derivative, with respect to x, gives 



■(l-«o(p)) 



a 2 Q,(u) 


/' dk 


- A;) - 6(1 


9x 2 


o V ? 


4 



<?'(!- p)-J(l-p) ... 6(l-k)6(l-p) 

- A n o{k) r 



8(g + 2fc x ) 



((Z + 2fc x ) 2 +4M 2 ' 



(B4) 



The derivatives only involve n^(k)(l — n„(p)), which re- 
sults in the square parenthesis. The third term in the in- 
tegrand is vanishing, since p = q = 1 implies q + 2k x = 
(note that p = k + q and q is along x) . Furthermore, 
we can change variable in the second term and cancel 
the no(p) contribution of the first term. Therefore, the 
square parenthesis simplifies to [S'(l — k) — 6(1 — k)]/A 
and the integration in dk is immediate. The final result 
coincides with the opposite of Eq. (|B3[) . Consider now 
the change of variable cos <f> — ^ k x , which gives 



d 2 Q q ( 



u 



dx 2 
f 1 8dk 



f 1 8k x dk x (q + 2k x ) 2 -Au 2 
-iq^T=kJ[(q + 2k x )*+4u 2 Y 

d (q + 2k x ) 2 - Au 2 



= 32 f ^^r^2( q + 2k 
J-i Q 



dk x {(q + 2k x ) 2 + 4u 2 ] 2 

(q + 2k x ) 2 - 12u 2 
\(q + 2k x ) 2 +4u 2 ] 3 ' 

(B5) 



0- 



where we integrated by parts in the second line. Notice 
that for q > 2 (see Appendix E]) j dk y n (k)(l — n (p)) = 



2^/1 — fc|. This establishes the equivalence of Eq. (|B5[) 
and (|Bip for this case. The equivalence of the two ex- 
pressions holds also at q < 2, as can be seen applying the 
change of variable k x — k x — q in the integration region 



-1 < k x < — | of (IB5p . Therefore, we conclude that 
d 2 Q q (u) d 2 Q q (u) _ d 2 Q q (u 



dg 2 



d\dg 



dx 2 



o 

(B6) 



This is in agreement with Eq. (|T6|) and leads to (|35[) . 
Appendix C: Physically alternate limit r s — > for 
finite g, X 

In this Appendix we analyze the r s — » limit of 
Eq. (|34|) for fixed values of g and %. In this situation 
the relevant integration region in dq is of order r s <C g, X 
around q — 0, and an 'anomalous' quadratic correction 
in g to the final result for the energy is obtained [see the 
discussion after Eq. (1331) ]. The calculation can be pat- 
terned after that carried out in the absence of spin-orbit 
coupling, as for instance done in Ref. In this case one 
can perform the integration in dq (up to an arbitrary up- 
per limit much larger than r s ) and extracts the coefficient 
of the r s In r s contribution by writing 



£ii(g,rs,x) 



r s In r s 
3V2(2tt) 4 J- 



[Q (u)] 3 du, (CI) 



which gives the standard result — 2 ^(10 — 37r)r s lnr s by 
using Qo(u) = AttR(u). Using for Qo(u) the expression 
of Eq. (|3"6]) instead, the result can still be obtained ana- 
lytically as in Ref. [33, i-e. 



£ii(g,r s ,x) 



In) 



6\/27r 



Ik 3 k 3 \ 



3k+k 2 



■F(k+,k- 



3/c — k j 

kJk\ 



F(k-,k + ) 



where we have defined the following function: 

V 2 % 2x 2 arccos - 
F(x, y) = A( x +y)-nx-4xE(l-^ 1 ) + ^= * , 

(C3) 

in terms of the elliptic function E(x), defined as in 
Ref. HH Notice that one can use the identity 

arccosh — 



/x-y 



fy—x 



when y > x. 
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